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General  Program  Overview 


The  objective  of  this  effort  was  to  enhance  sensitivity  and  develop  new  multi-mode 
formulations  for  radar  and  communications  in  the  presence  of  interference,  which  may  be 
coupled  across  the  dimensions  of  space,  slow-time,  fast-time,  and  polarization  according 
to  the  operating  mode.  In  particular,  the  effort  focused  on  waveform-diverse  emission 
schemes  such  as  various  physically  realizable  forms  of  MIMO  radar  and  polarization- 
diverse  radar,  adaptive  receive  processing  for  these  waveform-diverse  modes,  and  intra¬ 
pulse  radar-embedded  communications. 

Leveraging  the  notion  of  STAP  used  for  clutter  cancellation  in  airbome/space  radar, 
the  multiplicative  increase  in  degrees  of  freedom  afforded  by  multidimensional  coupling 
enables  expansion  to  other  dimensions  as  well.  Performing  this  coupling  on  transmit  also 
has  the  potential  to  improve  receive  detection/estimation  performance  via  the  resulting 
increased  diversity  of  the  transmit  signal.  This  effort  addressed  the  mathematical 
modeling  and  optimization  of  different  multidimensional  structures  for  the  realization  of 
new  physical  radar  emissions  which,  when  possible  based  on  available  test  equipment, 
were  also  evaluated  experimentally. 

New  adaptive  receiver  processing  schemes  were  also  likewise  developed  for  these 
multidimensional  emission  structures  which  were  also  evaluated  experimentally  based  on 
the  availability  of  measurements.  Finally,  this  effort  continued  the  theoretical 
development  of  a  scheme  for  waveform-diverse  embedding  of  low  probability  of 
intercept  communication  symbols  into  radar  clutter,  along  with  associated  symbol 
detection/decision  processing,  that  was  the  focus  of  a  previous  AFOSR  YIP  award.  It  is 
expected  that,  upon  transition,  these  waveform-diverse  schemes  will  serve  as  enabling 
technology  to  support  long-term  Air  Force  initiatives  to  address  an  ever  more  congested 
and  contested  RF  spectrum. 
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Section  I:  Multidimensional  Physical  Radar  Emissions 


As  outlined  in  [1],  which  resulted  from  collaborative  effort  under  the  NATO  SET-179 
and  SET- 182  research  task  groups  on  “Dynamic  Waveform  Diversity  &  Design”  and 
“Radar  Spectrum  Engineering  &  Management”,  respectively,  there  are  multiple 
challenges  facing  the  radar  community  as  a  result  of  an  increasingly  congested  and 
contested  RE  spectrum.  A  promising  approach  to  address  these  challenges  is  to  leverage 
advances  in  waveform  diversity  [2]  to  determine  how  to  better  exploit  the  multiple 
dimensions  in  which  a  radar  operates.  This  component  of  the  project  examined  how  to 
optimize  physical  radar  emissions  and  then  employed  this  same  framework  for  the  design 
of  physical  multidimensional  emissions. 

Optimization  of  Physical  Waveforms 

Within  the  context  of  the  polyphase-coded  EM  (PCEM)  framework  [3,4],  which 
established  the  means  to  convert  arbitrary  polyphase  codes  into  physically  realizable 
nonlinear  EM  (NEEM)  waveforms  amenable  to  the  distortion  induced  by  a  radar 
transmitter,  an  optimization  scheme  was  developed  to  search  this  very  high  dimensional 
space  [5].  Eor  a  code  with  N  +  I  chips,  where  each  could  take  on  one  of  M  different 
values  from  a  constant  amplitude  phase  constellation,  there  are  M  ^  possible  nonlinear 
EM  waveforms  that  could  be  generated  for  a  given  shaping  filter  g(t)  [4] .  Clearly  it  is  not 
possible  to  search  this  entire  space  if  either  M  or  N  are  moderately  large  (N  closely 
approximately  the  waveform  time-bandwidth  product  with  practical  values  from  roughly 
50  up  to  several  thousand).  Since  this  large  space  precludes  an  exhaustive  search  for  the 
global  optimum,  strategies  are  needed  (e.g.  [6,7])  to  determine  sufficiently  good  solutions 
according  to  the  desired  performance  measure  such  as  peak  sidelobe  level  (PSE)  or 
integrated  sidelobe  level  (ISE). 

Under  this  effort,  a  new  search  scheme  was  described  in  [5]  that  relies  on  the  fact  that 
PSE  and  ISE,  along  with  a  newly  defined  metric  denoted  as  frequency  template  error 
(ETE)  that  relies  on  the  Eourier  relationship  between  autocorrelation  and  power  spectral 
density  [5],  are  all  different  measures  of  the  waveform  autocorrelation  sidelobes.  Using 
this  relationship  it  is  possible  to  employ  a  greedy  search  strategy  that  avoids  the 
individual  local  minima  of  each  metric  by  alternating  between  them.  This  strategy, 
denoted  as  “performance  diversity”  is  a  twist  on  multi-objective  optimization  in  which, 
instead  of  seeking  to  balance  between  conflicting  metrics,  these  complementary  metrics 
provide  different  perspectives  on  the  same  underlying  criterion:  lower  range  sidelobes. 

Using  the  benchmark  of  the  PSE  bound  on  hyperbolic  EM  (HEM)  waveforms  [8] 
defined  as 


PSLhfm bound  =  -201ogio(5r)  -3  dB,  (1) 

which  is  a  function  of  time -bandwidth  product  (BT),  it  has  been  demonstrated  that 
optimized  PCEM  waveforms  can  surpass  this  benchmark  by  a  few  dB  for  a  given  N 
(which  well  approximates  BT).  It  will  be  demonstrated  later  that  new  ways  to  incorporate 
additional  dimensionality  into  the  waveform  without  increasing  BT  can  go  much  further. 

Eigures  1  and  2  illustrates  the  sidelobes  and  mainlobe  detail  for  use  of  the  individual 
metrics  (PSE,  ISE,  and  ETE)  and  the  performance  diversity  approach  along  with  the  EEM 
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waveform  for  comparison  when  BT  =  N  =6A  and  M  =  64.  In  this  case  the  PSL  benchmark 
from  (1)  is  found  to  be  -39.1  dB.  Table  1  provides  a  comparison  of  the  three  individual 
optimization  metrics,  performance  diversity,  and  LFM  in  terms  of  PSL,  ISL,  and  FTE  as 
well  as  mainlobe  resolution  defined  relative  to  the  baseline  LFM  waveform.  Note  that,  at 
a  cost  of  only  ~30%  broadening  of  range  resolution,  the  performance  diversity  optimized 
waveform  surpasses  the  HFM  bound  by  1.1  dB  (and  it  is  the  only  waveform  of  this  group 
to  even  get  close  to  doing  so). 

In  general  we  have  observed  that  such  optimized  FM  waveforms  that  can  surpass  the 
HFM  bound  largely  retain  the  chirp-like  structure  of  an  LFM  waveform  due  to  the 
consolidation  of  much  of  the  range-Doppler  ambiguity  into  the  well-known  ridge  that 
provides  Doppler  tolerance.  Such  waveforms  exhibit  the  generic  “sideways  S”  time- 
frequency  shape  that  is  typical  of  NLFM  waveforms  [8-12]  (and  summarized  in  [13]), 
albeit  with  seemingly  random  perturbations  that  arise  from  the  optimization  to  break  up 
the  coherency  that  otherwise  leads  to  higher  range  sidelobes.  In  fact,  some  of  these 
optimized  FM  waveforms  even  contain  short-term  wideband  components  that  are  quite 
similar  to  aspects  of  calls  generated  by  some  species  of  bats  [14],  which  unsurprisingly 
speaks  to  a  degree  of  optimality  in  the  bat  call. 


T 


Delay 

Figure  1:  Autocorrelation  of  the  optimized  FM  waveforms  using  the  PSL,  ISL,  and  FTE  metrics 
individually  and  the  performance  diversity  paradigm 
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Figure  2:  Autocorrelation  of  the  optimized  FM  waveforms  using  the  PSL,  ISL,  and  FTE  metrics 
individually  and  the  performance  diversity  paradigm  (mainlohe  close-up) 


Table  1:  Quantified  performance  for  optimized  FM  waveforms 


Optimization  Metrics 

LFM 

waveform 

PSL 

only 

!SL 

only 

FTE 

only 

Performance 

Diversity 

final  PSL  value  (dB) 

-13.5 

-21.1 

-27.0 

-31.3 

-4o.a 

final  ISL  value  (dB) 

-9.8 

-6.7 

-20.0 

-15.6 

-24.9 

final  FTE  value  (dB) 

-17.0 

-16.6 

-22.4 

-30.1 

-32.1 

relative  3  dB  resolution 

1.00 

i.m 

1.11 

1.26 

1.28 

relative  6  dB  resolution 

1.00 

1,08 

1.13 

1.30 

1.33 

It  is  also  described  in  [5]  how,  in  addition  to  designing  physical  waveforms  via 
optimization  of  the  underlying  coding  of  the  PCFM  framework,  one  may  extend  the 
optimization  to  encompass  the  impact  of  the  distortion-inducing  radar  transmitter  so  as  to 
ultimately  optimize  the  physical  emission  that  is  launched  into  free  space.  Figure  3 
depicts  this  arrangement  in  which  the  elements  of  the  discrete  code  are  collected  into  the 
vector  X,  the  operation  TpcpmI*}  converts  the  code  x  into  a  continuous  PCFM  waveform, 
and  rTx[*]  represents  the  distortion  induced  by  the  transmitter  onto  the  waveform  that 
subsequently  produces  the  free-space  emission.  It  is  shown  in  [5]  that  both  transmitter 
model-in-the-loop  (MiLo)  and  hardware-in-the-loop  (HiLo)  optimization  can  be 
performed.  The  latter  has  thus  far  been  demonstrated  with  low  power  systems  using  a 
Class  A  power  amplifier  (solid  state).  Future  work  will  examine  the  impact  of  high  power 
tube-based  amplifiers  (e.g.  klystron)  that  are  expected  to  induce  greater  distortion.  Such 
an  approach  may  provide  significant  benefit  to  ongoing  work  to  better  contain  the 
spectral  content  of  high-power  radar  emissions. 
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Figure  3:  Mathematical  representation  for  the  optimization  of  PCFM  radar  emissions 


Finally,  leveraging  previous  work  performed  in  collaboration  with  AFRL  [15] 
involving  the  optimization  of  a  non-repeating  form  of  nonlinear  FMCW  waveform  that  is 
optimized  on  a  segment-wise  basis,  the  notion  of  time-hopping  frequency  gaps  within  a 
radar  spectrum  was  explored  [16].  Because  it  is  essentially  a  form  of  FM  noise  radar,  as 
long  as  the  spectral  gaps  move  sufficiently  often  the  detrimental  increase  in  sidelobes  that 
would  otherwise  occur  is  largely  avoided.  Figure  4  illustrates  the  power  spectral  density 
for  individual  waveform  segments  within  which  fixed  spectral  gaps  reside  as  well  as  the 
aggregated  power  spectral  density  in  which  these  gaps  are  averaged  out  over  time.  The 
integrated  autocorrelation  over  a  coherent  processing  interval  (Fig.  5)  shows  that  moving 
the  spectral  gaps  more  often  leads  to  better  radar  sensitivity.  Future  work  will  explore 
how  this  approach  can  facilitate  the  incorporation  of  frequency-hopped  communications 
operating  in  tandem  with  the  frequency-hopped  radar  spectral  gaps. 


(a)  (b) 


Figure  4:  Power  spectral  density  for  (a)  10  sequentially  hopped  spectral  gaps  and  (b)  the  aggregated 
power  spectral  density  for  the  10  sequentially  hopped  gaps  and  100  randomly  hopped  gaps 
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Spatial  Modulation  of  Physical  Waveforms 

The  coded  physical  structure  of  the  PCFM  waveform  makes  it  attractive  as  a 
framework  from  which  to  construct  multidimensional  emission  schemes.  One  such 
scheme  that  was  developed  under  this  effort  denoted  as  “spatial  modulation”  [17,18] 
incorporates  an  additional  form  of  fast-time  coding  that  is  implemented  in  such  a  way  as 
to  control  the  spatial  steering  of  the  antenna  array  mainbeam  during  the  pulse  width.  In  so 
doing,  the  active  radar  illumination  can  mimic  a  similar  passive  attribute  of  the  human 
eye  (or  any  animal  possessing  fovea)  known  as  fixational  eye  movement  that  has  been 
found  to  enhance  visual  acuity  and  discrimination.  Defined  generally  as  the  “waveform 
diverse  array”,  this  formulation  also  subsumes  the  notion  of  the  frequency  diverse  array 
(e.g.  [19-23]). 

Defining  s(t ;  x)  as  the  PCFM  waveform  generated  by  the  underlying  code  x,  the 
actual  waveform  launched  from  the  mth  antenna  element  in  a  uniform  linear  array  (ULA) 
can  be  expressed  as 


sjt,0c,^,x^) 


1 


sit;x)b'”(f. 


(2) 


where 


' 

ft 

N 

Y 

b(t;xJ  =  Qxp- 

-j 

J^(0* 

X^„^(r-(n-l)7;) 

_n=\ 

dT  +  JSq 

A 
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is  the  FM  signal  that  steers  the  beam  in  fast  time  and  (•)”*  in  (2)  arises  from  the 
Vandermonde  structure  across  the  ULA.  The  vector  contains  the  sequence  of 
electrical  angle  phase  shifts  S2,  ■■■  to  steer  the  spatial  beam  in  fast-time  and  Aq  is 
the  initial  electrical  angle. 

The  far-field  emission  from  an  M-element  antenna  array  that  results  from  the 
combination  of  waveform  and  spatial  modulation  can  be  expressed  as 


1  (M-l)/2 

^  m=-(M-l)/2 


(4) 


for  0^  the  center  spatial  direction  of  the  beam  (relative  to  boresight),  d  the  antenna 
element  spacing  and  A  the  wavelength  (assuming  the  narrowband  assumption  holds),  with 
X  and  Xj,  omitted  for  compact  representation.  From  (4)  it  can  be  observed  that  the 
waveform-modulated  signal  changes  as  a  function  of  spatial  angle,  meaning  that  this 
formulation  is  a  physically  realizable  form  of  MIMO  radar.  Further,  by  maintaining  a 
coherent  beam  that  is  steered  in  fast-time,  the  spatially-modulated  emission  is  less  prone 
to  degradation  from  array  mutual  coupling  effects  that  can  otherwise  arise  for  arbitrary 
MIMO  emissions  [24,25]  and  that  could  potentially  damage  the  radar  [26,27]  due  to  high 
voltage  standing  wave  ratio  (VSWR). 

Using  the  emission  structure  of  (4),  the  delay-angle  dependent  received  signal  that 
results  from  the  reflection  of  the  spatially  modulated  radar  illumination  is,  for  the  mth 
antenna  element, 


y,„  (i.0^)=\[x(i.0) .  gd.e.e^) 


d0  +  u{t) . 


(5) 


In  (5),  the  term  x{t,0)  is  the  complex  scattering  as  a  function  of  delay  and  spatial  angle, 
u{t)  is  additive  noise,  and  *  represents  convolution.  Across  the  M  antenna  elements,  the 
collection  of  the  received  signal  defined  by  (5)  provides  a  delay-angle  coupled  response 
from  the  illuminated  environment,  where  the  coupling  naturally  realizes  a  multiplicative 
increase  in  degrees  of  freedom  for  receive  processing  in  the  same  manner  as  space-time 
adaptive  processing  (STAP)  for  airborne/space-based  radar.  The  implications  of  delay- 
angle  coupling  for  adaptive  receive  processing  are  demonstrated  in  Section  II,  where  the 
non-uniform  distribution  of  transmit  energy  as  a  function  of  spatial  angle  via  g  [t,  0, 0(, ) 

must  be  compensated  so  as  not  to  induce  a  noise  enhancement  effect. 

It  is  also  instructive  to  examine  the  impact  of  this  coupled  emission  when  non- 
adaptive  processing  is  performed  that  consists  of  standard  beamforming  and  pulse 
compression  (according  to  the  waveform-modulated  emission  in  each  spatial  direction). 
Figure  6  illustrates  the  non-adaptive  processing  responses  when  there  are  five  targets  in 
an  ‘X’  formation  in  delay-angle  for  each  of  four  different  illumination  schemes:  Case  1) 
standard  transmit  beamforming  (no  spatial  modulation).  Case  2)  a  first  null  to  first  null 
linear  spatial  sweep.  Case  3)  a  second  null  to  second  null  linear  spatial  sweep,  and  Case 
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4)  a  first  null  to  first  null  half-cycle  sinusoidal  spatial  sweep.  The  waveform  for  all  cases 
is  an  LFM  with  BT  =  N  =  200  and  the  center  direction  is  -  0°  (boresight). 

Because  all  but  the  center  target  reside  near  the  mainbeam  spatial  nulls,  the  standard 
transmit  beamforming  (Case  1)  only  clearly  identifies  the  center  target.  In  contrast,  the 
three  spatial  modulation  illumination  schemes  provide  much  better  visibility  of  all  five 
targets,  though  with  some  reduced  SNR  that  is  due  to  spreading  the  transmit  energy  over 
a  larger  spatial  angle  and  degraded  range  resolution  due  to  less  bandwidth  being  available 
for  any  given  direction  in  the  spatial  sweep.  It  is  clear  that  spatial  modulation  introduces  a 
new  tool  for  sensing  in  which  SNR  and  range  resolution  can  be  traded  for  enhanced 
target  visibility.  Examination  of  the  delay-angle  ambiguity  function  [17,18]  for  these  four 
emissions  (Fig.  7)  reveals  that  the  three  spatial  modulation  sehemes  also  naturally 
provide  enhanced  spatial  resolution  as  part  of  the  trade.  It  is  likewise  interesting  to  note 
the  impact  of  non-coherent  combining  of  multiple  pulses  to  exploit  the  phase  diversity  of 
different  targets  (Fig.  8),  where  the  target  configuration  is  now  clearly  visible  for  the 
three  spatial  modulation  emission  schemes. 
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Figure  6:  Five-target  ‘X’  scenario  for  Case  1)  standard  beamforming,  Case  2)  null-to-null  linear 
spatial  modulation,  Case  3)  double  null-to-null  linear  spatial  modulation,  Case  4)  null-to-null  balf- 

cycle  sinusoidal  spatial  modulation 
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Figure  7:  Peak-normalized  r  =  0 ,  0  =  0°  cut  of  angle-delay  ambiguity  function  (mainbeam  close-up) 
for  Case  1)  standard  beamforming,  Case  2)  null-to-null  linear  spatial  modulation,  Case  3)  double 
null-to-null  linear  spatial  modulation.  Case  4)  null-to-null  half-cycle  sinusoidal  spatial  modulation 
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Figure  8:  Four-pulse  non-coherent  integration  of  five-target  ‘X’  scenario  for  Case  1)  standard 
beamforming.  Case  2)  null-to-null  linear  spatial  modnlation.  Case  3)  double  null-to-null  Unear  spatial 
modnlation.  Case  4)  null-to-nnll  half-cycle  sinnsoidal  spatial  modulation 
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To  better  mimie  the  eye,  a  two-dimensional  (2D)  generalization  of  spatial  modulation 
was  demonstrated  in  [28]  for  a  planar  array.  The  spatially  modulation  FM  signal  from  (3) 
now  beeomes 


^.(^;x^,s)  =  exp 


VO 


n=\ 


dr  +  A 


x,0 


(6) 


and 


b(t;x  )  =  exp 


Vo 


n=\ 


dT  +  A 


z,0 


(7) 


for  X  and  z  the  azimuth  and  elevation  designations,  respeetively.  As  a  result,  (2) 
generalizes  to 


.s )  =  ^  ;  Xw )  [K  ;  Xx,s ))""  , ) )” 


(8) 


at  the  m  =  {m^,  m  )  antenna  element  so  that  the  far- field  emission  is  now 


(9) 


for  and  the  assoeiated  wave  numbers. 

For  eomparison,  eonsider  the  time-averaged  beampattems  (over  the  pulse)  shown  in 
Fig.  9  that  is  eomprised  of  standard  beamforming,  an  optimized  phase -dithered 
beamforming  that  is  eonstant  over  the  pulse,  and  spatial  modulation  that  traees  out  one 
eireular  rotation  during  the  pulse  with  a  radius  of  the  first  null  relative  to  the  eenter 
direetion.  For  two  elosely  spaeed  targets  residing  at  the  same  range.  Fig.  10  illustrates  the 
enhaneed  diserimination  eapability,  where  only  the  2D  spatial  modulation  ean  adequately 
diseriminate  the  two  targets  from  the  sidelobes  that  are  otherwise  indueed  by  either  form 
of  statie  beamforming.  As  with  the  ID  arrangement,  SNR  loss  and  range  resolution 
degradation  still  oeeur  with  the  2D  ease.  However,  it  is  not  diffieult  to  imagine  that,  from 
a  eognitive  sensing  standpoint,  sueh  trade-offs  may  be  useful  in  some  situations  in  mueh 
the  same  way  that  the  eye  performs  autonomie  adjustments  based  on  lighting  eonditions 
and  the  degree  of  aetive  attention  [29-31]. 
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Figure  9:  Aggregate  beampattern  in  Cartesian  coordinates  for  (a)  standard  beamforming;  (b)  phase- 
only  beamforming;  and  (c)  circular  spatial  modulation  with  center  direction 


Azimuth  0^  (degrees) 


Figure  10:  Received  response  for  (a)  standard  beamforming  (b)  phase-only  beamforming  and  (c) 
circular  spatial  modulation  with  targets  located  at  (0°,+ll.5°)  and  (+1.5°, -3°) 


Fast-Time  Polarization  Modulation 

Along  with  the  increased  dimensionality  of  fast-time  spatial  modulation  discussed 
above,  this  effort  also  explored  the  modulation  of  polarization  in  fast-time  via  a  similar 
mechanism  that  is  controlled  by  the  waveform.  In  so  doing,  the  waveform  and 
polarization  modulations  are  likewise  coupled,  thereby  increasing  the  dimensionality  of 
the  radar  emission,  with  applications  such  as  weather  radar  [32]  and  polarimetric 
synthetic  aperture  radar  (SAR)  [33]. 

The  simplest  way  in  which  fast-time  polarization  modulation  can  be  achieved  is  well- 
known  as  it  simply  involves  simultaneously  emitting  a  different  waveform  on  the 
horizontal  and  vertical  elements  of  a  dual-polarized  antenna  [34,35].  Since  the  desire  is 
often  to  separate  the  co-polarized  and  cross-polarized  components  of  the  received  signals 
on  the  horizontal  and  vertical  elements,  these  waveforms  should  have  low  cross¬ 
correlation,  with  the  most  common  waveforms  meeting  this  criterion  being  an  LFM  up- 
chirp  and  down-chirp.  In  this  case  the  instantaneous  polarization  during  the  pulse  will  lie 
somewhere  on  the  great  circle  shown  on  the  Poincare  sphere  in  Fig.  11. 
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Figure  11:  Feasible  polarizations  on  the  Poincare  sphere  for  two  waveforms  directly  connected  to 
horizontal  and  vertical  antenna  elements.  RCP:  right-hand  circular,  LCP:  left-hand  circular,  LHP: 

linear  horizontal,  LVP:  linear  vertical,  L-45:  linear  with  -45°  tilt,  L+45:  linear  with  +45°  tilt 

We  have  recently  shown  that  the  PCFM  waveform  structure  can  be  combined  with  a 
180°  hybrid  coupler  to  achieve  linear  amplification  using  nonlinear  components  (LINC) 
out  of  the  sum  channel  of  the  hybrid  as  a  means  to  slow  down  the  rapid  pulse  rise/fall¬ 
time  that  dominates  the  spectral  roll-off  of  the  radar  emission  [36].  In  this  case  an 
additional  form  of  coding  is  introduced  that  is  used  to  control  the  relative  phase  between 
two  otherwise  identical  waveforms  such  that,  when  the  two  waveforms  are  combined  in 
phase  within  the  hybrid  the  amplitude  can  be  controlled  (in  this  case  after  the  high  power 
amplifier  that  is  operated  in  saturation). 

This  same  idea  was  subsequently  applied  to  polarization  modulation  (Fig.  12),  where 
E  is  the  sum  channel,  A  is  the  difference  channel,  and  k  is  an  additional  phase  term  that 
dictates  the  particular  great  circle  that  can  be  traversed  on  the  Poincare  sphere.  The 
waveforms  5j(0  and  ^2(0  coded  such  that,  when  combined  in  the  hybrid,  they 
produce  the  resulting  ^^(0  and  Sy(t)  waveforms  that  are  emitted  from  the  horizontal 
and  vertical  antenna  elements,  respectively,  to  produce  the  desired  joint 
waveform/polarization  modulation  [37,38]. 


sXt)—  1  I  - 


s^(r) —  2  A  ■  --SvCO 


Figure  12:  180°  hybrid  for  polarization  control 


Mathematically,  the  result  from  the  hybrid  is 


(10) 


13 
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(11) 


1 

>^V  (0  =  ^  (.^2  (0  -  >^1  (0)  exp(7/t-) . 

For  the  examples  of  k  =  0°  and  k  =  -90°,  the  associated  great  circles  on  the  Poincare 
sphere  are  shown  in  Fig.  13.  Via  control  of  the  two  driving  waveforms  and  phase  term  k 
any  polarization  state  on  the  surface  of  the  Poincare  sphere  can  be  attained. 


RCP  RCP 


(a)  (b) 

Figure  13:  Feasible  polarizations  on  the  Poincare  sphere  using  the  180°  hybrid  for  a)  k  =  0°,  and 
b)  K  =  -90°.  RCP:  right-hand  circular,  LCP:  left-hand  circular,  LHP:  linear  horizontal,  LVP:  linear 
vertical,  L-45:  linear  with  -45°  tilt,  L+45:  linear  with  +45°  tilt 


The  waveforms  above  can  be  coded  to  enable  control  over  the  joint 
waveform/polarization  modulation  by  again  leveraging  the  PCFM  structure  as 


5i(0  =  exp); 


\g(Ty- 


vo 


N 

n=i 


dr  +  ^Q+y/Q 


(12) 


S2(t)  =  exp\j 


}g(r)  = 


VO 


N 

n=i 


dT  +  (/>Q-y/Q 


(13) 


where  «„  is  the  original  PCFM  waveform  coding,  serves  as  a  fast-time  polarization 
coding  about  the  particular  great  circle  of  the  Poincare  sphere  determined  by  the  value  of 
K,  and  ^0  and  y/^  dictate  the  initial  waveform  phase  and  polarization  state,  respectively. 

This  formulation  will  be  continued  to  be  explored  to  determine  the  utility  of  joint 
waveform/polarization  modulation  with  regard  to  various  sensing  applications,  metrics  by 
which  this  emission  scheme  may  best  be  evaluated  and  optimized,  and  for  broader 
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applicability  such  as  to  encode  communications  into  polarization  state.  In  Section  II  a 
new  adaptive  processing  scheme  is  described  that  is  applicable  to  dual-polarized 
emissions. 


Expanding  Waveform  Degrees  of  Freedom  beyond  the  Time -Bandwidth  Product 

The  time -bandwidth  product  (57)  of  a  waveform  is  a  measure  of  the  dimensionality 
of  the  waveform.  As  defined  in  (1),  the  peak  sidelobe  level  (PSL)  for  the  HFM  waveform 
has  a  known  bound  that  serves  as  a  benchmark  for  other  FM  waveforms.  In  this  context 
the  time  aspect  corresponds  to  the  pulsewidth  and  the  bandwidth  corresponds  to  the  3  dB 
bandwidth  (which  is  associated  with  matched  filter  range  resolution).  This  component  of 
the  project  focused  on  how  to  enhance  waveform  sensitivity  (reduce  sidelobes)  further 
beyond  the  HFM  bound  while  maintaining  the  same  57. 

The  PCFM  waveform  structure  described  above  enables  conversion  of  an  arbitrary 
polyphase  code  into  an  FM  waveform.  The  constellation  of  phase  values  from  which  a 
code  can  be  drawn  resides  on  the  unit  circle,  where  PCFM  converts  the  phase  change 
between  successive  code  values  into  a  continuous  phase  signal.  The  notion  of  “over¬ 
coding”  described  in  [39,40]  expands  this  arrangement  to  permit  1)  freedom  for 
additional  phase  trajectory  changes  during  a  chip  interval  and  2)  a  greater  amount  of 
phase  change  during  a  fixed  interval  as  long  as  the  aggregate  power  spectral  density  of 
the  waveform  does  not  exceed  a  specified  spectral  mask.  The  former  enables  even  finer 
“phase  dithering”  during  the  pulse,  an  attribute  that  has  been  observed  to  reduce  sidelobes 
by  breaking  up  their  coherence.  The  latter  enables  the  instantaneous  frequency  content  of 
the  waveform  to  increase  as  long  as  it  occurs  rarely  and  thus  conforms  to  the  desired 
power  spectral  density. 

Figure  14  illustrates  the  autocorrelation  of  three  optimized  PCFM  waveforms,  two  of 
which  employ  over-coding.  The  variable  L  defines  the  number  of  phase  “sub -transitions” 
during  a  chip  width  and  M  defines  the  amount  of  phase  transition  allowed  during  a  chip 
width  (+M7t:  so  M  >  1  exceeds  the  previous  polyphase-code -based  structure  defined  on 
the  unit  circle).  Compared  to  a  PSL  of  ^1.4  dB  without  over-coding  (L=l,  M=\),  which 
itself  surpasses  the  HFM  PSL  bound  of  -39.1  dB  by  2.3  dB,  the  (L=8,  M=l)  and  (L=8, 
M=2)  over-coded  waveforms  achieve  PSL  values  of  -44.1  dB  and  -52.0  dB,  respectively, 
with  the  latter  necessitating  a  spectral  mask  to  prevent  the  3  dB  bandwidth  from 
increasing.  The  sensitivity  enhancement  arises  from  increased  design  degrees  of  freedom 
that  are  evident  from  the  now  more  rapid  variation  of  instantaneous  frequency  (see  Fig. 
15). 

The  power  spectral  density  is  shown  in  Fig.  16  (with  the  (7=8,  M=\)  case  replaced  by 
an  optimized  (7=4,  M=2)  case  to  show  a  different  combination  of  parameters).  In 
particular,  note  that  the  M  =  2  cases  exhibit  some  spectral  spreading  (or  “spectral  fuzz”) 
in  the  roll-off  region,  though  the  imposed  Gaussian  power  spectral  density  mask  is  still 
maintained.  Thus  this  notion  of  over-coding  raises  the  fundamental  question  of,  for  a 
prescribed  spectral  mask  that  dictates  roll-off  and  mainlobe  width  (to  define  range 
resolution),  just  how  low  can  the  range  sidelobes  be  driven  for  a  constant  amplitude 
waveform?  Note  that  our  group  recently  showed  [41]  that  a  small  amount  of  amplitude 
modulation  (with  associated  mismatch  loss  of  ~0.25  dB)  can,  if  jointly  optimized  with  the 
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FM  waveform,  achieve  a  physical  emission  with  sidelobes  below  -108  dB  in  simulation 
and  -83  dB  demonstrated  experimentally. 


Delay 

Figure  14:  Autocorrelation  comparison  for  over-coded  optimized  waveforms  (L=8,  M=1  and  L=8, 
M=2)  relative  to  the  non-over-coded  optimized  waveform  (L=l,  M=l)  for  BF  s  64 


Figure  15:  Normalized  instantaneous  frequency  for  the  three  waveforms 
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Figure  16:  Power  spectral  density  for  over-coded  waveforms  -  note  the  (L  =  4,  M  =  2)  waveform  has 
replaced  the  (L  =  8,  M  =  1)  waveform  here  to  show  a  different  combination 


Alternative  to  over-coding,  in  [42,43]  a  higher-order  framework  for  PCFM 
waveforms  was  described.  Defining  the  instantaneous  waveform  phase  as  (jiit),  a 
combined  first,  second,  and  third  order  phase  formulation  can  be  expressed  as 


t  T  t'  r  t  t 

^3(0  =  |||c(r")dr"dr'dr-i-||Z7(r')  dr' dr  +^a{T)  (14) 

000  00  0 


where 


a(r)  = 

n=\ 


(15) 


is  the  first-order  coded  function  according  to  the  values  in  the  sequence  and  initial 
phase  ^3  0 ,  which  collectively  represent  the  original  PCFM  framework.  The  term 


Kr’)=Z^„g2(r’-(n-l)r^) 


«=1 


(16) 
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is  the  second-order  coded  function  representing  time-varying  chirp-rate  according  to  the 
values  in  the  sequence  while 


n=\ 


(17) 


is  the  third-order  coded  function  representing  time-varying  chirp  acceleration  according 
to  the  values  in  the  sequence  c„ .  The  terms  gi(r),  g2(^')’  associated 

shaping  filters  for  each  order.  Thus  (14)  involves  three  different  coding  structures  that 
can  be  manipulated  to  define  a  given  FM  waveform.  It  has  been  found  that  the  first  and 
second  order  terms  yield  the  greatest  benefit  to  designing  good  radar  waveforms.  A 
smaller  additional  benefit  is  obtained  from  the  third  order  term  with  improvement 
continuing  to  diminish  rapidly  as  the  order  increases  further.  For  example,  Fig.  17  depicts 
the  autocorrelations  of  waveforms  that  were  jointly  optimized  across  the  different  orders. 
When  the  and  F'  orders  were  jointly  optimized  a  PSL  value  of  ^8.4  dB  was  attained. 
When  the  3“^  order  component  was  incorporated  for  joint  optimization  the  PSL  only 
improved  to  -48.7  dB  (a  difference  of  0.3  dB). 


Delay 


Figure  17:  Autocorrelations  of  jointly  optimized  waveforms 


It  is  interesting  to  consider  the  instantaneous  frequency  for  these  two  waveforms  as 
shown  in  Fig.  18.  The  familiar  NLFM  structure  is  apparent,  though  a  close  examination 
reveals  the  small  dithering  that  is  typical  of  PCFM  to  break  up  sidelobe  coherence.  Future 
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work  will  investigate  the  combination  of  this  higher-order  framework  with  the  above 
over-coding  scheme. 
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Section  II:  Multidimensional  Adaptive  Receive  Processing 


In  Section  I  multiple  new  physical  emission  schemes  were  described.  While  the 
straightforward  matched  filter  may  be  applied  to  the  resulting  received  echoes,  it  is  useful 
for  many  applications  to  perform  some  form  of  optimal  or  adaptive  receive  processing. 
However,  previous  methods  for  optimal  [44]  and  adaptive  pulse  compression  [45]  were 
based  on  a  (non-physical)  polyphase  code  framework  that  is  not  directly  extensible  to 
physical  waveforms. 

This  component  of  the  project  focused  on  the  generalization  of  the  adaptive  pulse 
compression  (APC)  algorithm  to  permit  application  to  FM  waveforms  which  constitute 
the  majority  of  high-power  radar  emissions  being  used  in  practice.  The  framework  of  the 
problem  can  be  stated  as  follows: 

1.  The  phase  function  for  an  FM  waveform  is  continuously  changing  and  thus 
the  waveform  must  be  sampled  at  a  rate  high  enough  to  provide  acceptable 
fidelity  (greater  than  the  3  dB  bandwidth  to  account  for  sufficient  spectral  roll¬ 
off  since  a  time-limited  pulse  has  theoretically  infinite  bandwidth). 

2.  Such  “over-sampling”  with  respect  to  3  dB  bandwidth  results  in  higher 
dimensionality  that 

a.  increases  computational  cost, 

b.  leads  to  ill-conditioning  effects  when  sidelobes  are  suppressed  (for 
optimal  and  adaptive  processing), 

c.  induces  a  range  super-resolution  condition  that  produces  spurious 
sidelobe  peaks. 

3.  Addressing  the  above  issues  then  enables  application  of  adaptive  processing  to 
physical  emissions  that  possess  spatial  diversity,  polarization  diversity,  etc. 
via  incorporation/compensation  of  the  salient  physical  features. 

In  [46,47]  the  APC  structure  has  been  generalized  to  account  for  these  physical 
effects  through  the  combination  of  1)  an  MVDR-based  structure  to  prevent  over- nulling 
[48],  2)  a  polyphase  decomposed  representation  of  multiple  delay-shifted  versions  of  the 
over-sampled  waveform  [49],  and  3)  a  “beam  spoiling”  modification  to  prevent  the 
degradation  associated  with  super-resolution  [50]  (relative  to  the  matched  filter).  Figure 
19  shows  experimentally  measured  results  using  an  LFM  waveform  taken  from  the  roof 
of  Nichols  Hall  on  the  University  of  Kansas  campus.  These  receive-processed  results 
include  the  new  version  of  APC  along  with  the  traditional  matched  filter  (MF)  and  two 
new  instantiations  of  Least-Squares  optimal  mismatch  filtering  (MMF)  suited  to  FM 
waveforms.  The  MMF  results  clearly  outperform  the  MF,  though  the  APC  results  yield 
the  most  significant  enhancement.  Using  an  optimized  FM  waveform  from  Section  I 
reveals  similar  performance  (Fig.  20),  though  the  difference  between  the  different  receive 
processing  schemes  is  less  significant  due  to  use  of  a  far  better  waveform. 
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Figure  19;  Pulse  compressed  response  using  LFM  waveform  (simultaneous  transmit  and  receive 
using  separate  antennas  so  direct  path  dominates) 
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Figure  20:  Pulse  compressed  response  using  optimized  FM  waveform  (simultaneous  transmit  and 
receive  using  separate  antennas  so  direct  path  dominates) 
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The  version  of  APC  amenable  to  physical  FM  waveforms  has  also  been  extended  to 
incorporate  physical  MIMO  emissions.  Subsuming  previous  work  on  direction  of  arrival 
estimation  [51,52]  and  a  prior  instantiation  developed  for  (non-physical)  polyphase  codes 
[53],  this  approach  is  denoted  as  space-range  adaptive  processing  (SRAP)  as  it  jointly 
operates  in  the  range  and  spatial  angle  domains  [54].  For  the  spatial  modulation  scheme 
described  in  Section  I,  Fig.  21  depicts  the  matched  filter  and  SRAP  responses, 
respectively,  when  five  proximate  targets  are  illuminated  by  a  spatially  modulated 
emission  that  linearly  sweeps  in  both  frequency  and  angle.  It  is  observed  that  SRAP 
suppresses  delay-angle  sidelobes  and  enhances  resolution  in  both  the  range  and  spatial 
domains.  This  enhanced  discrimination  capability  tends  to  incur  a  cost  in  terms  of  SNR 
so  such  an  approach  is  likely  to  be  most  applicable  when  SNR  is  sufficiently  high  to 
permit  enhanced  discrimination.  The  mathematical  derivation  and  implementation  of 
SRAP  can  be  found  in  [54]. 
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Figure  21:  Five  targets  illuminated  by  spatially  modulated  emission  using  (a)  delay-angle  matched 

filtering  and  (b)  the  SRAP  algorithm 


This  waveform-based  adaptive  framework  has  also  been  developed  for  application  to 
simultaneous  dual-polarized  radar  emissions  as  a  means  to  disambiguate  the  co-polarized 
and  cross-polarized  reflections  at  the  radar  receiver.  The  polarimetric  APC  (PAPC) 
algorithm  derived  in  [37,38]  realizes  a  different  adaptive  range-domain  filter  for  each 
individual  range  cell  and  received  polarization  state.  An  experimental  demonstration  of 
this  capability  is  illustrated  in  Fig.  22.  An  up-chirp  and  down-chirp  having  BT  =  100  were 
concurrently  emitted  from  the  horizontal  and  vertical  polarized  elements  using  a  dual-pol 
Vivaldi  antenna.  An  identical  antenna  was  used  on  receive  to  enable  simultaneous 
transmit  and  receive.  As  with  the  above  results,  the  direct  path  becomes  the  dominant 
signal  that  masks  many  of  the  smaller  proximate  echoes.  A  dihedral  corner  reflector  was 
present  about  97  m  from  the  radar  testbed  which  shows  up  clearly  in  the  cross-polarized 
responses. 

The  blue  trace  in  the  data  is  the  matched  filter  response  that  demonstrates  the  usual 
cross-correlation  one  expects  from  using  two  waveforms  and  which  masks  the  smaller 
scatterers.  The  ultra-low  sidelobe  (ULS)  response  [41]  (green  trace)  is  obtained  from 
using  an  optimized  waveform  that  is  emitted  from  the  H  and  V  antennas  at  separate  times 
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to  obtain  an  estimate  of  ground  truth.  We  find  that  the  PAPC  algorithm  operating  on  the 
echoes  from  the  up/down-chirp  waveforms  (red  trace)  exhibits  performance  that  is  nearly 
identical  to  the  time-separated  optimized  ULS  response  (e.g.  the  small  echo  from  Nichols 
Hall  is  now  observable).  Thus  one  can  conclude  that  PAPC  provides  an  adaptive  solution 
to  enable  simultaneous  dual-polarized  radar  emissions. 


0  20  40  60  80  100  120 


Figure  22:  Experimental  demonstration  of  adaptive  separation  of  radar  echoes  generated  by 
simultaneous  dual-polarized  radar  emissions,  (a)  HH  response;  (b)  HV  response;  (c)VH  response;  (d) 
VV  response.  Matched  Filter  (blue),  PAPC  (red)  and  ULS  (green). 
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Section  III:  Intra-Pulse  Radar-Embedded  Communications 


This  final  aspect  of  the  project  continued  the  investigation  of  a  form  of  covert 
communication  embedded  via  a  tag/transponder  into  radar  backscatter  (i.e.  clutter)  that 
was  the  main  focus  of  a  previous  AFOSR  YIP  project  [55-57].  Recently  in  [58]  the 
practical  aspects  of  this  problem  were  delineated  to  provide  a  framework  for  continued 
development,  both  theoretical  and  experimental  (see  Table  2). 

Table  2:  Practical  design  constraints  for  intra-pulse  radar-embedded  communications 

1 .  Radar,  tag,  and  desired  receiver  are  not  synchronized  (no  control  channel 
available). 

2.  Communication  symbols  must  be  determined  independently  at  the  tag  and 
desired  receiver  based  on  the  observed  radar  illumination  (must  address 
possible  timing  uncertainty  and  multipath  differences). 

3.  The  surrounding  clutter  structure  within  which  the  symbols  are  to  be 

covertly  embedded  is  not  known  exactly,  though  the  clutter  response  can  be 
viewed  as  a  random  signal  convolved  with  the  radar  waveform. _ 

4.  Symbols  should  be  designed  to  be  partially  correlated  with  the  clutter 

response  (in  the  range  domain)  but  still  separable  from  one  another  on 
receive. _ 

5.  Symbols  should  have  a  temporal  and  spectral  footprint  commensurate  with 
the  radar  waveform  to  ensure  proper  masking  by  the  clutter  and  thus  remain 
low  probability  of  intercept  (LPI). 

6.  Symbol  design  must  occur  in  a  reasonable  time  frame  (while  radar 
illumination  persists). 


It  is  analytically  demonstrated  in  [58]  that  spectral  shaping  of  the  set  of 
communication  symbols  (Fig.  23)  according  to  the  spectral  footprint  of  the  radar 
emission  (and  thus  resulting  clutter  response)  enhances  the  desired  low  intercept 
probability  (LPI)  attribute  of  the  symbols  while  retaining  good  communication 
performance  in  terms  of  bit  error  rate.  Specifically,  the  processing  gain  at  a  desired 
receiver  that  has  knowledge  of  the  possible  symbols  is  evaluated  for  different  symbol 
design  strategies  (Fig.  24).  These  same  symbol  strategies  were  likewise  evaluated 
according  to  the  processing  gain  that  would  be  obtained  at  a  hypothetical  “worst  case” 
intercept  receiver  that  possesses  partial  clairvoyant  knowledge  of  the  symbols  (i.e.  their 
time  width  and  bandwidth). 

The  comparison  of  these  two  processing  gains  then  reveals  the  gain  advantage  of  the 
desired  receiver  thereby  realizing  the  underlying  design  metric  that  captures  the 
fundamental  goodness  of  a  given  symbol  design  approach  (e.g.  Fig.  25).  Of  the  various 
symbol  design  strategies  considered,  it  was  the  Shaped  Water  Filling  (SWF)  approach 
that  was  found  to  provide  the  best  all-around  performance  in  terms  of  minimizing 
communication  error  rate  and  ensuring  LPI  [58]. 
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Digital  Freq. 

Figure  23:  Spectral  content  for  radar-embedded  communication  symbols  where  R+N:  radar 
emission  with  noise,  DP:  dominant  projection  symbols,  SDP:  shaped  dominant  projection  symbols, 
SWF:  shaped  water  fiUing  symbols,  and  DSSS:  direct -sequence  spread  spectrum  symbols  (for 
baseline  comparison). 


Figure  24:  Comparison  of  processing  gain  as  a  function  of  symbol  dimensionality  m  for  30  dB  radar 
clutter  to  noise  ratio  using  “best”  observed  receiver  filtering  (DF:  decorrelating  filter  or  LDF:  loaded 
decorrelating  filter)  for  each  symbol  set. 


DISTRIBUTION  A:  Distribution  approved  for  pubiic  reiease 


25 


Figure  25:  Probability  of  detection  (solid)  and  intercept  (dashed)  for  different  dimensionalities  of 
SWF  symbols.  The  SNR  difference  for  the  same  dimensionaUty  demonstrates  the  gain  advantage. 
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